Importance of Correlation Effects on Magnetic Anisotropy in Fe and Ni 
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We calculate magnetic anisotropy energy of Fe and Ni by taking into account the effects of strong 
electronic correlations, spin-orbit coupling, and non-collinearity of intra-atomic magnetization. The 
LDA+U method is used and its equivalence to dynamical mean-field theory in the static limit is 
emphasized. Both experimental magnitude of MAE and direction of magnetization are predicted 
correctly near U = 1.9 eV, J = 1.2 eV for Ni and U = 1.2 eV, J = 0.8 eV for Fe. Correlations 
modify the one-electron spectra which are now in better agreement with experiments. 



The calculation of the magneto-crystalline anisotropy 
energy (MAE) 0-H] of magnetic materials containing 
transition-metal elements from first principles electronic 
structure calculations is a long-standing problem. The 
MAE is denned as the difference of total energies with 
the orientations of magnetization pointing in different, 
e.g., (001) and (111), crystalline axis. The difference is 
not zero because of spin-orbit effect, which couples the 
magnetization to the lattice, and determines the direc- 
tion of magnetization, called the easy axis. 

Being a ground state property, the MAE should 
be accessible in principle via density functional theory 
(DFT) ||7| . Despite the primary difficulty related to the 
smallness of MAE (~ 1 /LteV/atom), great efforts to com- 
pute the quantity with advanced total energy methods 
based on local density approximation (LDA) combined 
with the development of faster computers, have seen suc- 
cess in predicting its correct orders of magnitudes [|| |l2[ . 
However, the correct easy axis of Ni has not been pre- 
dicted by this method and the fundamental problem of 
understanding MAE is still open. 

A great amount of work has been done to understand 
what is the difficulty in predicting the correct axis for Ni. 
Parameters within the LDA calculation have been varied 
to capture physical effects which might not be correctly 
described. These include (i) scaling spin-orbit coupling 
in order to enlarge its effect on the MAE |j,[l(J, (ii) cal- 
culating torque to avoid comparing large numbers of en- 
ergy ]irj|j , (iii) studying the effects of the second Hunds 
rule in the orbital polarization theory pdj |, (iv) analyz- 
ing possible changes in the position of the Fermi level by 
changing the number of valence electrons JT3], (v) us- 
ing the state tracking method p3| , and (iv) real space 
approach p4J . 

In this paper we take a new view that the correla- 
tion effects within the d shell are important for the mag- 
netic anisotropy of 3d transition metals like Ni. These 
effects are not captured by the LDA but are described 
by Hubbard-like interactions presented in these systems 
and need to be treated by first principles methods jlq] . 

Another effect which has not been investigated in the 
context of magnetic anisotropy calculations is the non- 
collincar nature of intra-atomic magnetization |lq] . It is 



expected to be important when spin-orbit coupling and 
correlation effects come into play together. In this arti- 
cle we show that when we include these new ingredients 
into the calculation we solve the long-standing problem 
of predicting the correct easy axis of Ni. 

We believe that the physics of transition metal com- 
pounds is intermediate between atomic limit where the 
localized d electrons are treated in the real space and 
fully itinerant limit where the electrons are described 
by band theory in k space. A many-body method in- 
corporating these two important limits is the dynamical 
mean-field theory (DMFT) 0. The DMFT approach 
has been extensively used to study model Hamiltonian of 
correlated electron systems in the weak, strong and inter- 
mediate coupling regimes. It has been very successful in 
describing the physics of realistic systems such as tran- 
sition metal oxides and, therefore, is expected to treat 
properly the materials with d or / electrons. 

Electron-electron correlation matrix U lll2l3 ~ u = 
(mi?7i3 \vc\ mim/±) <5 SlS2 <5 S3S4 for d orbitals is the quantity 
which takes strong correlations into account. This matrix 
can be expressed via Slater integrals F^ l > , i = 0, 2, 4, 6 in 
the standard manner. The inclusion of this interaction 
generates self-energy E 7l72 (io> n , k) on top of the one 
electron spectra. Within DMFT it is approximated by 
momentum independent self-energy E 7l72 (iw n ). 

A central quantity of the dynamical mean field theory 
is the one-electron on-site Green function 
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where H® (k) is the one-electron Hamiltonian stan- 
dardly treatable within the LDA. Since the latter already 
includes the electron-electron interactions in some aver- 
aged way, we subtract the double counting term Vd c P8| . 
The use of realistic localized orbital representation such 
as linear muffin-tin orbitals p9[ | leads us to include over- 
lap matrix 7l72 (&) into the calculation. 

The DMFT reduces the problem to solving effec- 
tive impurity model where the correlated d orbitals are 
treated as an impurity level hybridized with the bath of 



conduction electrons. The role of hybridization is played 
by the so-called bath Green function defined as follows: 
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Solving this impurity model gives access to the self- 
energy £ 7l72 (iw„) for the correlated electrons. The one- 
electron Green function (|lj) is now modified with new 
S 7l72 (iu! n ), which generates a new bath Green function. 
Therefore, the whole problem requires self-consistency. 

In this paper we confine ourselves to zero temperature 
and make an additional assumption on solving the im- 
purity model using the Hartree-Fock approximation. In 
this approximation the self-energy reduces to 
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is the average occupation matrix for the cor- 



related orbitals. The off-diagonal elements of the occu- 
pancy matrix are not zero when spin-orbit coupling is 
included |2C|j . The latter can be implemented following 
the prescription of Andersen p9[ or more recent one by 
Pederson J2l]] . 

In the Hartree-Fock limit the self-energy is frequency 
independent and real. The self-consistency condition of 
DMFT can be expressed in terms of the average occu- 
pation matrix: Having started from some n 7l72 we find 
S 7l72 according to (ra). Fortunately, the computation of 
the on-site Green function (pi) needs not to be performed. 
Since the self-energy is real, the new occupancies can 
be calculated from the eigenvectors of the one-electron 
Hamiltonians with £ 7l72 — Vdc added to its dd block. The 
latter can be viewed as an orbital-dependent potential 
which has been introduced by the LDA+U method |15|j . 

The LDA+U method has been very successful com- 
pared with experiments at zero temperature in ordered 
compounds. By establishing its equivalence to the static 
limit of the DMFT we see clearly that dynamical mean- 
field theory is a way of improving upon it, which is crucial 
for finite temperature properties. 

In this work we study the effect of the Slater param- 
eters -Fo, -F2 and F4 on the magnetic anisotropy energy. 
Slater integrals can be linked to intra-atomic repulsion 
U and exchange J obtained from LSDA supercell pro- 
cedures via U = F° and J = (F 2 + F 4 )/14. The ratio 
F 2 /F A is to a good accuracy a constant ~ 0.625 for d 
electrons |22J| . The MAE is calculated by taking the 
difference of two total energies with different directions 
of magnetization (MAE=E(111) - £(001)). The total 
energies are obtained via fully self consistent solutions. 
Since the total energy calculation requires high precision, 
full potential LMTO method pj| has been employed. 
For the k space integration, we follow the analysis given 
by Trygg and co-workers [[fl[ and use the special point 
method [E4J with a Gaussian broadening pa| of 15 mRy. 
The validity and convergence of this procedure has been 
tested in their work [fill. For convergence of the total 



energies within desired accuracy, about 15000 fc-points 
are needed. We used 28000 fc-points to reduce possi- 
ble numerical noise, where the convergency is tested up 
to 84000fc-points. Our calculations include non-spherical 
terms of the charge density and potential both within 
the atomic spheres and in the interstitial region |23|] . All 
low-lying semi-core states are treated together with the 
valence states in a common Hamiltonian matrix in order 
to avoid unnecessary uncertainties. These calculations 
are spin polarized and assume the existence of long-range 
magnetic order. Spin-orbit coupling is implemented ac- 
cording to the suggestions by Andersen fL9|| . We also 
treat magnetization as a general vector field, which re- 
alizes non-collinear intra-atomic nature of this quantity. 
Such general magnetization scheme has been recently dis- 
cussed pq . 

To incorporate the effects of intraatomic correlations 
on the magnetocrystalline anisotropy energy, we have to 
take into account the intra-atomic repulsion U and the 
intraatomic exchange J. It is important to perform the 
calculations for fixed values of magnetic moments which 
themselves show a dependency on U and J. Since the 
pure LSDA result (C/=0, and J=0) reproduces the ex- 
perimental values for magnetic moments in both Fe and 
Ni fairly well, we have scanned the U — J parameter space 
and have obtained the path of U and J values which hold 
the theoretical moment constant, following the approach 
ofRef. |||. 

We now discuss our calculated MAE. We first test our 
method in case of LDA (U = J = 0). To compare with 
previous calculations, we turn off the non-collinearity of 
magnetization which makes it collinear with the quanti- 
zation axis. The calculation gives correct orders of mag- 
nitude for both fee Ni and bec Fe but with the wrong 
easy axis for Ni, which is the same result as the previous 
result |ll[] . Turning on the non-collinearity results in a a 
larger value of the absolute value of the MAE (2.9 /J,eV) 
for Ni but the easy axis predicted to be (001) which is 
still wrong. The magnitude of the experimental MAE of 
Ni is 2.8 fieV aligned along (111) direction |27|. 

We now describe the effect of correlations, which is 
crucial in predicting the correct axis of Ni (see Fig. fa). 
We walked along the path of parameters U and J which 
hold the magnetic moment to 0.6 /j,b- The MAE first 
increases to 60 fieV (U = 0.5 eV, J = 0.3 eV) and then 
decreases. While decreasing it makes a rather flat area 
from U = 1.4 eV, J = 0.9 eV to U = 1.7 eV, J = 1.1 eV 
where MAE is positive and around 10 /ieV. After the 
flat area, the MAE changes from the wrong easy axis to 
the correct easy axis. The correct magnetic anisotropy 
is predicted at U — 1.9 eV and J — 1.2 eV. The change 
from the wrong easy axis to the correct easy axis occurs 
over the range of SU ~ 0.2eU, which is of the order of 
spin-orbit coupling constant (~ O.leU). 

For Fe, the MAE is calculated along the path of U and 
J values which fixes the magnetic moment to 2.2 [is- 
At U = eV and J = eV, the MAE is 0.5 fieV. It 
increases as we move along the contour in the direction of 



increasing U and J. The correct MAE with the correct 
direction of magnetic moment is predicted at U = 1.2 eV 
and J = 0.8 el/. 

Notice that the trends in the values of U and J that are 
necessary to reproduce the correct magnetic anisotropy 
energy within LDA + U are similar to the values used 
to describe photoemission spectra of these materials pq] 
within DMFT. The values of the parameters U and J, 
are basis-set dependent, and method dependent, but the 
values of U used in our LDA+ U calculation are within 
leV of those used in Q. Since DMFT contains the 
graphs which screen the on-site interaction which are 
ommitted in the LDA+ U functional, a larger value of U 
is needed to produce the correct moment in DMFT. 

We find direct correlation between the dependency of 
the MAE as a function of U-J and the difference of mag- 
netic moments (Am = — (m(lll) — m(001)) behaving 
similarly (see Fig. [l]). For Ni the difference increases till 
U = 0.4 eV and J = 0.2 eV, then decreases. While de- 
creasing it makes a flat area from U = 0.9 eV and J = 0.6 
to U = 1.7 eV and J = 1.1 eV. After the flat area, 
the difference decreases rapidly. For Fe, the difference 
of magnetic moments slightly fluctuates till U = 0.7 eV 
and J = 0.5 eV and then decreases till U = 1.0 eV and 
J = 0.7 eV . 

This concurrent change of MAE and the difference of 
magnetic moments suggests why some previous attempts 
based on force theorem |lj] failed in predicting the cor- 
rect easy axes. Force theorem replaces the difference of 
the total energies by the difference of one-electron ener- 
gies. In this approach, the contribution from the slight 
difference in magnetic moments does not appear and, 
therefore, is not counted in properly. Unfortunately, we 
could not find any experimental data of magnetic mo- 
ments with different orientations to the desired precision 
(10~ 4 /is) to compare with. 

We now present implications of our results on the cal- 
culated electronic structure for the case of Ni. One im- 
portant feature which emerges from the calculation is 
the absence of the X2 pocket (see Fig. ||). This has 
been predicted by LDA but has not been found exper- 
imentally | p9| . The band corresponding to the pocket 
is pushed down just below the Fermi level. This is ex- 
pected since correlation effects are more important for 
slower electrons and the velocity near the pocket is rather 
small. It turns out that the whole band is submerged un- 
der the Fermi level. We also find that the removal of the 
X2 point is near the point U = 1.9 eV and J = 1.2 eV. 
For comparison, the corresponding band is just above the 
Fermi level at U = 1.9 eV and J = 1.1 eV forming a tiny 
pocket. This strengthens the connection between MAE 
and the absence of X2 pocket. 

There has been some suspicions that the incorrect po- 
sition of the X2 band within LDA was responsible for the 
incorrect prediction of the easy axis within this theory. 
Daalderop and coworkers |l2| removed the X 2 pocket by 
increasing the number of valence electrons and found the 
correct easy axis. We therefore conclude that the absence 



of the pocket is one of the central elements in determining 
the magnetic anisotropy, and there is no need for any ad- 
hoc adjustment within a theory which takes into account 
the correlations. 

We now describe the effects originated from (near) de- 
generate states close to the Fermi surface. These have 
been of primary interest in past analytic studies poLpTJ . 
We will call such states degenerate Fermi surface crossing 
(DFSC) states. The contribution to MAE by non-DFSC 
states comes from the fourth order perturbation. Hence 
it is of the order of A 4 . The energy splitting between 
DFSC states due to spin-orbit coupling is of the order 
of A because the contribution comes from the first order 
perturbation. Using linear approximation of the disper- 
sion relation e(fcA), the relevant volume in fc-space was 
found of the order A 3 . Thus, these DFSC states make 
contribution of the order of A 4 . Moreover, there may be 
accidentally DFSC states appearing along a line on the 
Fermi surface, rather than at a point. We have found 
this case in our LDA calculation for Ni. Therefore the 
contribution of DFSC states is as important as the bulk 
non-DFSC states though the degeneracies occur only in 
small portion of the Brillouin zone. 

The importance of the DFSC states leads us to com- 
parative analysis of the LDA and LDA+U band struc- 
tures near the Fermi level. In LDA, five bands are cross- 
ing the Fermi level at nearly the same points along the 
TX direction. Two of the five bands are degenerate for 
the residual symmetry and the other three bands acciden- 
tally cross the Fermi surface at nearly the same points. 
There are two sp bands with spin up and spin down, 
respectively. The other three bands are dominated by 
d orbitals. In LDA+U, one of the d bands is pushed 
down below the Fermi surface. The other four bands 
are divided into two degenerate pieces at the Fermi level 
(see Fig. g): Two symmetry related degenerate d± bands 
and two near degenerate sp-j- and spi bands. In LDA, 
we found that two bands are accidentally near degener- 
ate along the line on the Fermi surface within the plane 
TXL. One band is dominated by d± orbitals. The other 
is dominated by d± orbitals near X and by s± orbitals 
off X. In LDA+U, these accidental DFSC states disap- 
pear (see Fig. |j). 

As we have seen, the on-site repulsion U reduces the 
number of DFSC states along TX direction. Based on 
the tight -binding model, the importance of DFSC states 
has been shown. We see that strong correlations reduce 
number of DFSC states in TX direction and remove the 
near degenerate states on TXL plane. We conclude that 
the change of DFSC states is another important element 
that determines the easy axis of Ni. 

To conclude, we have demonstrated that it is possible 
to perform highly precise calculation of the total energy 
in order to obtain both the correct easy axes and the 
magnitudes of MAE for Fe and Ni. This has been ac- 
complished by including the strong correlation effects via 
taking intra-atomic repulsion and exchange into account 



and incorporating the non-collinear magnetization. In 
both Fe and Ni, both U and J take physically accept- 
able values consistent with the values known from atomic 
physics. The calculations performed are state of the art 
in what can currently be achieved for realistic treatments 
of correlated solids. Further studies should be devoted 
to improving the quality of the solution of the impurity 
model within DMFT and extending the calculation to 
finite temperatures. 
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FIG. 1. The magnetocrystalline anisotropy 

energy MAE = £'(111) - £(001) (circle) and the difference 
of magnetic moment Am = m(001) — m(lll) (square) for 
Ni (top) and for Fe (bottom) as functions of U. The experi- 
mental MAEs are marked by arrows for Fe (1.4 jieV) and Ni 
(—2.8 neV). The values of exchange parameter J for every 
value of U are chosen to hold the magnetic moment of 0.6 /is 
in Ni and 2.2 \xb in Fe 



FIG. 2. Calculated Fermi Surface of Ni with the correla- 
tion effects taken into account. The solid and dotted lines 
correspond to majority and minority dominant spin carriers. 
Dominant orbital characters are expressed. Both experimen- 
tally confirmed Y5 pocket and L neck can be seen. The X2 
pocket is missing, which is in agreement with experiments. 
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